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ABSTRACT 

We present the analysis of 23 published rotation curves of disk galaxies belonging 
to the Ursa Major group of galaxies, with kinematics free of irregularities. The rotation 
curves are analysed in the context of MOND (Modified Newtonian Dynamics). We add 
an extra component to the rotation curve fits, in addition to the stellar and gaseous 
disks: a speculative halo of constant density made of, e.g., neutrinos, which would 
solve the bulk of the problem currently faced by MOND on rich galaxy clusters scales. 
We find that this additional unseen mass density is poorly constrained (as expected 
a priori, given that a neutrino halo never dominates the kinematics), but we also find 
that the best-fit value is non-zero: p = 3.8 x 10~^^ g cm~'^, and that a zero-density 
is marginally excluded with 87% confidence; also, the 95% confidence upper limit for 
the density is p = 9.6 x 10~^^ g cm""^. These limits are slightly above the expectations 
from the Tremaine-Gunn phase space constraints on ordinary 2 eV neutrinos, but in 
accordance with the maximum density expected for one or two species of 5 eV sterile 
neutrinos. 

Key words: galaxies: kinematics and dynamics - dark matter - galaxies: spiral - 
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1 INTRODUCTION 

The current dominant paradigm in cosmology is that galac- 
tic dark matter is made of non-baryonic weakly-interacting 
massive particles, the so-called cold dark matter (CDM). 
However, on galaxy scales, observations are at variance with 
a sizeable list of CDM predictions (e.g., Klypin et al. 1999; 
Gentile et al. 2004, 2005, 2007a, 2007b, 2007c; Simon et al. 
2003; Famaey & Binney 2005; Kuzio de Naray et al. 2006), 
whilst one observes that the gravitational field is closely re- 
lated to the baryonic distribution in spiral galaxies (e.g. 
McGaugh et al. 2007; Famaey et al. 2007a) through the 
modified Newtonian dynamics (MOND) relation of Milgrom 
(1983): in MOND, for gravitational accelerations a below 
ao 10~*cms~^ the effective a approaches (gNOio)^^^ where 
gTN is the usual Newtonian gravitational field. This leads to 
remarkable fits of galactic rotation curves over five decades 
in mass ranging from tiny dwarfs (e.g.. Gentile et al. 2007a, 
2007c) through early-type spirals (e.g., Sanders & Noorder- 
meer 2007) to massive ellipticals (Milgrom & Sanders 2003) , 
without resorting to galactic dark matter. MOND also ac- 
counts very well for observed scaling relations of galaxies 
(e.g., Sanders & McGaugh 2002; McGaugh 2004, 2005). 



Moreover, the mysteriously small but non-zero cosmologi- 
cal constant also finds a natural answer in some co-variant 
versions of MOND (e.g., Zhao 2007). 

However, despite all these successes, MOND fails to 
completely explain away the mass discrepancy in clusters of 
galaxies (Aguirre, Schaye & Quataert 2001; Sanders 2003, 
2007; Pointecouteau & Silk 2005; Angus, Famaey & Buote 
2007a): consequently, these objects still require substantial 
amounts of non- luminous matter in MOND. The recent 
weak-lensing observations of the "bullet cluster" (Clowe et 
al. 2006, Angus et al. 2007b) have moreover constrained 
this non-luminous component to be coUisionless, e.g. dark 
baryons in some non-dissipative form (MACHOs or small 
dense clumps of cold gas). Another possibility, suggested by 
Sanders (2003), is that ordinary neutrinos have a mass of 
the order of 2 eV (near the present experimental limit) , and 
thus contribute significantly to the mass budget of galaxy 
clusters. These hot dark matter particles are fermionic and 
have a maximum phase space density set by equilibrium in 
the early Universe. Reexpressing this phase space limit from 
Tremaine & Gunn (1979), and assuming 3 species of ordi- 
nary neutrinos of nearly equal mass, Sanders (2003) noted 
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that in a relaxed cluster of virialised gas temperature T, if 
the neutrino fluid has the same temperature as the gas, the 
maximum neutrino density is: 

3/2 



g cm 



(1) 



where rrii, is the ordinary neutrino mass. 

It has indeed been shown that galaxy clusters scaling 
relations, including the luminosity-temperature relation, are 
then naturally reproduced in MOND (Sanders 2007), while 
such a non-baryonic component of moderately heavy neutri- 
nos helps fitting the angular power spectrum of the Cosmic 
Microwave Background in relativistic MOND theories (e.g., 
Skordis et al. 2006). However, it has recently been shown 
that such ordinary 2 eV neutrinos would fail to explain 
the mass discrepancy in X-ray emitting groups of galaxies, 
where the X-ray temperature falls below 2 keV (Angus et 
al. 2007a), because the maximum neutrino density would 
be too low. This notion might be related to the fact that, 
due to their unrelaxed nature, the temperature of hot gas 
in groups is not representative of the actual temperature 
of the neutrino fluid. Alternatively, this could rather hint 
at the presence of extra baryonic mass in the form of cold 
dense clumps of cold gas (e.g. Pfenniger, Combes & Mar- 
tinet 1994; Pfenniger & Combes 1994), or at the existence 
of one or two species of light sterile neutrinos with masses of 
a few eV (e.g. Maltoni et al. 2004, Maltoni & Schwetz 2007). 

In any case, if neutrinos do contribute significantly to 
the mass budget of rich galaxy clusters, it is then likely that, 
in the context of MOND, individual galaxies and galaxy 
groups possess extensive low-density neutrino halos. These 
were suggested by Sanders (2007) to possibly affect the weak 
lensing properties of galaxies on scales of hundreds of kpc. In 
this Letter, we rather investigate the possible contribution 
of neutrinos on the rotation curves of individual galaxies re- 
siding in groups in the context of MOND. To achieve this, 
we reanalyze the rotation curves in the Ursa Major galaxy 
group (Sanders & Verheijen 1998) by adding as a new pa- 
rameter a constant neutrino density. 



2 "MOND + NEUTRINOS" MODELS 

In MOND, the transition from the Newtonian regime (where 
a = (/n) to the deep MOND regime (where a — (^nho)^''^) is 
regulated by an interpolation function /x(a;), where x — a/ao 
and a = g-^^/ ^,{x). fJ.{x) should have the following asymptotic 
behaviours: j-i{x) 1 for x » 1 and ij.{x) — > a/ao for x « 
1. Traditionally the following functional form of fi{x) was 
used: ^j.{x) — ^2-^1/2 1 but recently the preferred form has 
become the "simple" n{x) proposed by Famaey & Binney 
(2005): 



fi{x) 



1 + x 



(2) 



This form of fi{x) has been shown to give disk galax- 
ies' rotation curves fits which are as good as the standard 
/i(a;), with realistic mass-to-light [M / L) ratios (Famaey et 
al. 2007a). Also, contrary to the standard ^(x), it gives rea- 
sonable disk and bulge M /L ratios in early-type disk galax- 
ies (Sanders & Noordermeer 2007). 

To motivate a 2 eV neutrino component in galaxies and 
galaxy groups, we note that massive neutrinos are hot dark 



matter (decoupled from electrons and the radiation at rela- 
tivistic speed), but that these thermal relics have presently 
cooled to very low speed. Indeed, since their momentum 
comparable to that of CMB photons today, their speed is 
expected to be ~ 30 km s"^ x (4/ll)~^/^ = 21 km s~\ 
The factor (4/11)"^'''^ arises from entropy conservation ar- 
guments in the early Universe (e.g. Weinberg 1972). They 
can however be heated up to 100-1000 km s~^ speed by 
falling into the gravitational potentials of group and clus- 
ters of galaxies. 

Eq. [1] could be generalised to situations where the 
X-ray temperature is unknown or unapplicable. Adopting 



1 keV(- 



from Sarazin (1986), we can convert 



^ 400 km a-l • 

gas temperature to the 1-dimensional equilibrium velocity 
dispersion a. The Tremaine-Gunn limit for three species of 
ordinary neutrinos (and their antiparticles) of nearly equal 
mass can then be recasted as: 



4.8 X 10-27 g cm- 



1/3 



\2 ev) 



-4/3 



400 km s~^ 

rp .1/2 



(3) 



.1 keV 

Note that 4.8 x 10"^'' g cm"^ = 7 x 10"^ Mq pc"^. 
3 DATA 

The Ursa Major group is an ideal example to perform this 
kind of rotation curve analysis. Galaxy groups are typically 
dominated by spiral galaxies for which the density profiles 
can be measured by analysing their rotation curves. Fur- 
thermore, such groups are usually difi^use and of relatively 
low mass. Consequently, the external field effect in galaxy 
groups is less important than in more massive galaxy clus- 
ters (Bekenstein & Milgrom 1984, Famaey, Bruneton & Zhao 
2007b, Angus & McGaugh 2007, Wu et al. 2007). The Ursa 
Major galaxies do not form a rich, "classical" virialised clus- 
ter, such as the ones for which 2 eV neutrinos have been in- 
voked (Sanders 2003, 2007) as a solution to the MOND mass 
discrepancy problem. It is known for a while that such non 
X-ray emitting groups do not exhibit a mass discrepancy 
problem in MOND (e.g., Milgrom 2002), but within the un- 
certainties due notably to the unrelaxed nature of the group, 
a low density extra mass component cannot be excluded. 

Note that an individual spiral galaxy is expected to 
be surrounded by a very shallow and extended neutrino 
halo (Sanders 2007), according to the phase-space contraint 
which corresponds to a temperature below 0.1 keV (i.e. 
typically a ^ 120 km for an individual galaxy, see e.g. 
Battaglia et al. 2005 for the Milky Way). Such halos, sur- 
rounding individual galaxies in a group, would merge with 
the halo of the group as a whole, for which the equivalent 
"virialised gas temperature" is T=0.14 keV, corresponding 
to a dispersion of 150 km . However, given that the Ursa 
Major group is not virialised (Sanders & Verheijen 1998), 
this value can only be taken as a very rough estimate. 

Hereafter, we model the neutrino component around 
each galaxy as the halo of the whole Ursa Major group, 
with an isotropic constant density. Given the uncertainties 
(the actual neutrino temperature, their density profile, the 
physical distances of the galaxies from the group centre). 
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it is a realistic assumption since a halo of ordinary 2 eV 
neutrinos is expected to feature a constant density core of 
several hundreds kpc (Sanders 2003) . More detailed neutrino 
density distributions will be the subject of future studies. We 
then just fit rotation curves the usual way in MOND, but 
adding this extra component. 

One of the great advantages of the galaxies of this sam- 
ple is that they can be assumed to lie at approximately the 
same distance, since they belong to the same group. Con- 
trary to Sanders & Verheijen (1998), who assumed a dis- 
tance of 15.5 Mpc, here we took the distance to be 18.6 
Mpc, the value calculated by TuUy & Pierce (2000) from 
the TuUy-Fisher relation and Cepheid-determined distances 
to nearby galaxies; we also note that Sakai et al. (2000) de- 
rive a slightly higher value of the distance, 20.7±3.2 Mpc, 
consistent with 18.6 Mpc. 

The data that we used in the present paper were taken 
from Sanders & Verheijen (1998). The authors present HI 
rotation curves of 30 galaxies belonging to the Ursa Major 
group. The derivation of the rotation curves is presented 
in Verheijen (1997) and Verheijen & Sancisi (2001); they 
were derived from HI observations performed at the WSRT 
(Westerbork Synthesis Radio Telescope). The original sam- 
ple included 62 galaxies brighter than Mb ~ —16.5, but 
32 of them were not considered by Sanders & Verheijen for 
their kinematical analysis because they either had a too low 
inclination, they were lacking HI, they were interacting with 
neighbouring galaxies, or had too few independent points in 
the rotation curve. Their final sample of 30 galaxies spans 
a large range of Hubble types. In the present analysis, how- 
ever, we only consider the 23 galaxies which are free of kine- 
matical disturbances. The stellar mass distribution is de- 
rived from NIK (Tf'-band) photometry (TuUy et al. 1996), 
and the linear resolution of the rotation curve is 0.75 kpc. 



4 RESULTS AND DISCUSSION 

We attempted to constrain the neutrino density, using the 
MOND hypothesis, by fitting the 23 rotation curves dis- 
cussed in the previous section. 

The rotation curves were fitted the usual way in MOND, 
but adding an extra component in the fits due to the pu- 
tative neutrino halo. We used the simple fi{x) function (see 
Section [21). The global parameters of the fit were ao and the 
additional (constant) density p^: they both were constrained 
to be the same for all the galaxies of the sample. We note 
that keeping ao as a free parameter can compensate any 
errors on the (fixed) distance. Then, the stellar (i^'-band) 
mass-to- light (M/L) ratio was left as an additional free pa- 
rameter, separately for each galaxy. The uncertainties on the 
fitted parameters were derived from the statistics of the 
global fit performed on the 23 rotation curves of the sample. 

We find a 1-a allowed range for between 0.8 x lO"^'^ 
g cm~^ and 7.4 x 10~^^ g cm~^, the best-fit value being = 
3.8 X 10"^^ g cm~"^. The corresponding best-fit value of ao 
is Q.72±0.09 X 10"* cm s"^. In Fig.[T]we plot the confidence 
levels in parameter space: as shown by Begeman, Broeils & 
Sanders (1991), fits with variable ao and fixed distance D 
are essentially identical to fits with fixed ao and variable D, 
because the product aoD^ comes in the computation of the 
MOND circular velocity. Note that higher distances suggest 



a lower neutrino density. Also, a negative neutrino density 
cannot be measured, thus any fit of the density leads to a 
non-negative value. 

As expected a priori, the constraints on are very 
weak, as the neutrino halo never dominates the kinematics 
due to its very low density. Fig. [5] shows the overall effect 
of adding a neutrino component: at large radii, where the 
neutrinos contribution is larger, the residuals of the fits be- 
come smaller than without neutrinos. In the outer parts, 
both model rotation curves are below the observations, but 
the curves with neutrinos are closer to the observations. The 
wave-like behaviour of the residuals in Fig. [2] will be studied 
in a future paper: different p-functions might explain some 
systematic residual of the rotation curves (Sanders & Ver- 
heijen 1998; Milgrom & Sanders 2007). However, using the 
standard p-function does not get rid of the wave-like be- 
haviour of Fig. (2] We also note that p^ = is excluded with 
87% confidence. 

Let us now check whether the neutrino fits are consis- 
tent with the Tremaine-Gunn limit and the global observa- 
tion of Ursa Major. Plugging the 1-a range of our rotation 
curve- determined value (0.8 x 10"^^ to 7.4 x 10"^'^) in the 
Tremaine-Gunn limit for ordinary 2 eV neutrinos, we predict 
that a = 220-460 km s"\ or equivalently that T = 0.3-1.3 
keV. These should be the typical equilibrium dispersion and 
temperature for characterizing the overall potential which 
the 2 eV neutrinos are in. This equilibrium velocity disper- 
sion is probably marginally consistent with the global poten- 
tial of the Ursa Major group: no X-ray emission has been 
detected, but it is observed to have a galaxy velocity disper- 
sion of aobs = 148 ± 6 km s"^ (TuUy 1987). However, Ursa 
Major is unrelaxed, so (Tobs is likely an under-estimate of 
the potential. On the other hand, the lack of X-ray implies 
the intergalactic gas in Ursa is most certainly below 1 keV. 
Adopting a = 148 km s^^, the data is slightly pushing for 
a somewhat higher than 2 eV neutrino mass, which could 
be reminiscent of the fact that ordinary 2 eV neutrinos are 
not found to explain away the discrepancy in X-ray emitting 
groups (Angus et al. 2007a). 

Let us finally note that a possible effect on all these 
results is the external field effect (EFE), discussed in Section 
(2] galaxies in clusters are embedded in an external field much 
higher than the typical value of O.Olao in which the Milky 
Way is embedded (Famaey et al. 2007b, Wu et al. 2007), so 
even when their internal acceleration is much smaller than 
ao they are actually not in a MOND regime. However, as 
argued in Sanders & Verheijen (1998), the EFE is expected 
to be of order 0.02ao in the Ursa Major group, meaning that 
we can neglect it as a first order approximation. 



5 CONCLUSIONS 

We analysed 23 rotation curves of disk galaxies with regu- 
lar kinematics belonging to the Ursa Major group of galax- 
ies, presented by Sanders & Verheijen (1998). The rotation 
curves are fitted within the MOND (Modified Dynamics, 
Milgrom 1983) context, where traditionally only the visible 
matter (gas and stars) is needed to fit rotation curves. How- 
ever, given that MOND shows a mass discrepancy in rich 
galaxy clusters which can be mainly solved by adding a halo 
of 2 eV ordinary neutrinos (Sanders 2003), we add a third 
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component to the MOND fits, with the same constant den- 
sity for all galaxies of the sample. It is however important 
to bear in mind that this extra mass component in our fits 
does not necessarily have to be ordinary neutrinos. 

The fits give a rather unconstrained best-fit value of the 
additional constant density: = 3.8I3 Q x g cm"''. 

Such weak constraints are to be expected since a neutrino 
halo is never the dominant component in the kinematics. 
However, a zero additional density is marginally excluded, 
with 87% confidence. These results show that neutrino (or 
any other dark matter) effects are worth being modelled in 
MOND in the future, with better data and larger samples. 

We however comment that there is a slight tension be- 
tween the rotation curve-determined neutrino density and 
the Tremaine-Gunn limit for three species of ordinary 2 eV 
neutrinos, assuming that the velocity dispersion is a good 
estimator of the neutrino fluid temperature (which is un- 
sure in such an unrelaxed cluster). Indeed, from Eq.(l), and 
for T = 0.14 keV, the maximum neutrino density for UMa 
should be 2.5 x 10"^** (m,./2eV)* g cm~^. For ordinary 2 eV 
neutrinos, this density is a factor 15 smaller than our best-fit 
value. A similar tension is observed in fitting weak-lensing 
data of high-z X-ray clusters (Takahashi & Chiba 2007, An- 
gus et al. 2007b), temperature profiles of X-ray emitting 
groups (Angus et al. 2007a), and velocity dispersions of 
galaxies at the center of rich clusters (Richtler et al. 2007). 
This means that the narrow window of ordinary neutrino 
mass of 2-2.2 eV is presently tightly squeezed by opposite 
constraints from astronomical data and particle physics ex- 
periments. Actually, this tension could rather hint at the 
existence of one, two or three species of light sterile neutri- 
nos (Maltoni & Schwetz 2007). For three species of sterile 
neutrinos with a 5 eV mass (and ordinary neutrinos with a 
negligible mass), a mass rriv = 5 eV can be directly plugged 
in Eq. (1). In the case of two species of sterile neutrinos, 
the maximum density must then be divided by a factor 1.5, 
and in the case of one species by a factor 3. This means 
that for two species of 5 eV sterile neutrinos, the maximum 
density for a temperature of 0.14 keV is 6.5 x 10"'^^ g cm"'', 
and for one species 3.25 x lO"^'^ g cm~^, densities much 
closer to our rotation curve-determined values. These one 
or two species of 5 eV sterile neutrinos could additionally 
solve the problem of the MOND missing mass at the center 
of high-z X-ray clusters (Takahashi & Chiba 2007, Angus et 
al. 2007b), in X-ray emitting groups (Angus et al. 2007a), 
and in individual galaxies at the center of clusters (Richtler 
et al. 2007). 

Lastly, we note that doing similar rotation curve work 
in the much hotter Virgo and Coma clusters might reveal 
a stronger effect of the neutrinos, expected to be 100 times 
denser in Coma than in Ursa Major. It is thus of interest to 
fundamental physics to push for accurate resolved kinematic 
data of galaxies in the Virgo and Coma clusters. 
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Figure 1. 68% (solid curve) and 95% (dashed curve) confidence 
levels. D is the distance and the cross marks the best-fit value. 
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